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Abstract 

We investigate various classes of Gauge Mediated Supersymmetry Breaking models 
and show that the Next-to-Minimal Supersymmetric Standard Model can solve the 
^-problem in a phenomenologically acceptable way. These models include scenarios 
with singlet tadpole terms, which are phenomenologically viable, e.g., in the presence 
of a small Yukawa coupling < 10~^. Scenarios with suppressed trilinear j4-terms 
at the messenger scale lead naturally to light CP-odd scalars, which play the role of 
pseudo i?-axions. A wide range of parameters of such models satisfies LEP constraints, 
with CP-even Higgs scalars below 114 GeV decaying dominantly into a pair of CP-odd 
scalars. 
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1 Introduction 



The mediation of supersymmetry breaking to the observable sector via supersymmetric 
gauge interactions (GMSB) has aheady been proposed during the very early days of super- 
symmetric model building [1,2]. The essential ingredients of this class of models are a 
sequestered sector containing a spurion or a dynamical superfield X, whose F-component 
Fx does not vanish (there could exist several such fields). In addition, a messenger sector 
(fii exists, whose fields have a supersymmetric mass M, but a mass splitting between its 
scalar /pseudoscalar components due to its coupling to Fx- They carry Standard Model 
gauge quantum numbers such that the messengers couple to the Standard Model gauge su- 
permultiplets. Possible origins of supersymmetry breaking in the form of a nonvanishing 
Fx component can be O'Raifeartaigh-type models [2], models based on no-scale supergrav- 
ity [3,4] or Dynamical Supersymmetry Breaking [5-7]. 

If supersymmetric gauge interactions would be the only interactions that couple the 
visible sector with the messenger/sequestered sector, the phenomeno logically required fi 
and Bfi terms of the Minimal Supersymmetric Standard Model (MSSM) would be difficult 
to generate. The simplest solution to this problem is the introduction of a gauge singlet 
superfield S and a superpotential including the \SHuHd term, which has been used in early 
globally [8] and locally supersymmetric [9] models. 

Let us point out a possible connection between gravity mediated supersymmetry breaking 
and GMSB-like models [3, 4]: standard gravity mediated supersymmetry breaking within the 
MSSM requires Giudice-Masiero-like terms (depending on the Higgs doublets) in the Kahler 
potential [10] in order to generate the // and terms (see [11] for a possible 5-dimensional 
origin of such terms). Given a possible source for such terms, one can replace the Higgs 
doublets by the messengers of GMSB models and proceed as in the usual analysis of gauge 
mediation. The advantage of such models is that no other gravity mediated source of super- 
symmetry breaking as scalar or gaugino soft masses is required; such sources of supersym- 
metry breaking are frequently absent in higher dimensional setups. On the other hand, the 
solution of the standard /x-problem for the Higgs doublets still requires the introduction of a 
singlet S. Then one is also led to the scenario considered in this paper, the Next-to- Minimal 
Supersymmetric Standard Model (NMSSM) with gauge mediated supersymmetry breaking. 

In order to generate a sufficiently large vacuum expectation value of the scalar component 
S* of S" (and hence a sufficiently large effective /i term jji^ff = A (5*)), the singlet superfield S 
should possess additional Yukawa interactions with the messenger /sequestered sector. Then, 
an effective potential for S with the desired properties can be radiatively generated. 

Note that the so-called singlet tadpole problem [12] is absent once the original source of 
supersymmetry breaking is of the F-type [3, 4, 13]. On the contrary, singlet tadpole diagrams 
can now generate the desired structure of the singlet effective potential [3,4], triggering a 
VEV of S. If the singlet couples at lowest possible loop order to the messenger /sequestered 
sector such that tadpole diagrams are allowed, a mild version of the singlet tadpole problem 
reappears, since the coefficients of the corresponding terms linear in S are typically too large. 
This milder problem can be solved under the assumption that the involved Yukawa coupling 
is sufficiently small - however, it does not need to be smaller than the electron Yukawa 
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coupling of the Standard Model (see below). 

In the meantime, quite a large number of models involving GMSB and at least one gauge 
singlet, that generates an effective fi term, have been studied [14-19]. They differ in the 
particle content of the messenger /sequestered sector, and include sometimes more than one 
gauge singlet superfield. 

The purpose of the present paper is the investigation of a large class of models obtained af- 
ter integrating out the messenger /sequestered sector (including possibly heavy singlet fields). 
It is assumed that the remaining particle content with masses below the messenger scale M 
is the one of the NMSSM. 

The couplings and mass terms of the NMSSM are obtained under the following assump- 
tions: 

- no interactions between the Higgs doublets Hu, Ha and the messenger /sequestered sector 
exist apart from supersymmetric gauge interactions; then no MSSM-like fi ot Bfi terms are 
generated after integrating out the messenger /sequestered sector; 

- the gauge singlet superfield S has Yukawa interactions with the messenger /sequestered 
sector. As a result, various soft terms and S'-dependent terms in the superpotential can be 
generated after integrating out the messenger /sequestered sector. 

Under the only assumption that the original source of supersymmetry breaking is Fx 
and that the messengers have a mass of the order M > ^/Fx, superspace power counting 
rules allow to estimate the maximally possible order of magnitude of the generated masses 
and couplings. 

In general, these masses and couplings will comprise nearly all possibilities consistent 
with gauge invariance (see Section 2), leading to the general NMSSM. However, many of 
these mass terms and couplings can be much smaller, or absent, than indicated by the 
power counting rules (but never larger), if the corresponding diagrams involve high loop 
orders, small Yukawa couplings, or are forbidden by discrete or (approximate) continuous 
symmetries. 

In the next Section, we will parametrize the mass terms and couplings of the general 
NMSSM, and estimate their (maximally possible) radiatively generated order of magnitude 
with the help of superspace power counting rules. Section 3 is devoted to a phenomenological 
analysis of three different scenarios, which are defined by particular boundary conditions for 
the NMSSM parameters at the messenger scale, and Section 4 contains our conclusions. 

2 Results of superspace power counting rules 

The class of models investigated in this paper is defined by a superpotential 

W = XSH^Hd + ^S' + W{S,X,^i,. ..) + ... , (2.1) 

where W{S,X,ipi, . . .) denotes the couplings of 5* to the messenger/ sequestered sector, and 
we have omitted the standard Yukawa couplings of and Ha- No MSSM-like /z-term is 
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assumed to be present. Due to a coupling X(pi(pi in W , a non-vanishing Fx-component 

Fx = (2.2) 

induces a mass term 

^-^«+^;^) (2-3) 

which gives opposite contributions to the squared masses of the real and imaginary compo- 
nents of the scalar components of the messengers l^i. Since we assume no direct couplings 
of S to X, this constitutes the only original source of supersymmetry breaking. 

After integrating out the messenger /sequestered sector, the remaining effective action for 
the light superfields $ (the fields S*, Hu-, Hd, ... of the NMSSM) is necessarily of the form 

^c, /" d^eMD^,D^,$J,X,X) , (2.4) 

i 

where the relevant terms are obtained after the replacement of at least one superfield X 
by its F-component Fx- The maximally possible orders of magnitude of the coefficients q 
can be obtained by dimensional analysis: if a function /j is of a mass dimension [M]'^^, the 
corresponding coefficient q has a mass dimension [M]'^~'^f . As long as c// > 2 (which will be 
the case) , q will typically depend on the mass of the heaviest particle running in the loops to 
the appropriate power, and subsequently we identify this mass M with a unique messenger 
scale M,i^f,gg. 

We are aware of the fact that models exist where the Cj depend on several mass scales 
Mi] however, it is always trivial to identify a mass scale M such that Cj are bounded from 
above by M'^~'^f . Also, in the particular case df = 2, Cj can involve large logarithms; these 
depend on whether the VEV Fx is "hard" (i.e. generated at a scale A much larger than M) 
or "soft", i.e. generated by a potential involving terms of the order of M. In the first case, 
logarithms of the form ln(A^/M^) can appear in q. 

In the present situation (no interactions between the Higgs doublets Hu, Hd and the 
messenger/sequestered sector) possible supercovariant derivatives Da, Da inside ft do not 
lead to terms that would otherwise be absent; for this reason we will omit them in our 
analysis. (Here, we will not discuss the radiatively generated gaugino masses and scalar 
masses for the gauge non-singlets, but concentrate on the NMSSM specific effects.) To 
lowest loop order we can use the underlying assumption that only the singlet superfield 5* 
has direct couplings to the messenger /sequestered sector (however, see Fig. 1 below). The 
first terms that we will investigate are then of the form 

J2c^ [ d'eu(S,'s,X,X) . (2.5) 

Below we list all relevant terms with this structure. Given an expression of the form 
(12.51) . the generated S- and Fg-dependent terms can be obtained by the replacements 

X = M + ^W, S = S + e'^Fs . (2.6) 
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Due to the coupling X(pi(pi, the supersymmetry conserving mass M of the messengers (pi 
can be identified with the value of the scalar component of X. Loop factors like (IGtt^)"^ 
and model dependent Yukawa couplings are not explicitly given, but we indicate the powers 
of m (which follow from the powers of Fx) and M (which follow from dimensional analysis). 

The possible operators /j and the corresponding contributions to the scalar potential are 
then given by: 

SX + h.c. : m^Fs + h.c. (2.7) 

4 

SXX + h.c: —S + m'^Fs + h.c. (2.8) 

2 

^ 771 

SSX + h.c: —{SF* + h.c.) + FsF* (2.9) 

4 9 

777 777 

—SS* + —{SF* + h.c.)+FsF* (2.10) 

2 

^^T^ 777 

^5X + /i.c. : — (5F5 + /I.C.) (2.11) 

4 2 
^ ^ ^ ^ 777 777 

^5XX + /i.c. : + ^"^^^ + ^-c- (2-12) 

Operators with higher powers of X or X do not generate new expressions, and operators 

with higher powers of or generate negligible contributions with higher powers of M in 
the denominator (recall that we are assuming M >j m). 

The terms ~ FsFg in (12.91) and (12.101) only account for a correction to the wave func- 
tion normalization of the superfield S, which can be absorbed by a redefinition of S. The 
remaining terms can be written as an effective superpotential AW and additional contribu- 
tions AVsoft to the soft terms of the general NMSSM. To this end, the terms SFg + h.c. in 
(12.91) and (I2.10p have to be rewritten using the expression derived from the superpotential 

F*s = \H^Hd + kS'^ + ... (2.13) 

where the dots stand for terms of higher order in the loop expansion. We parametrize the 
effective superpotential AW and the soft terms AVsoft of the general NMSSM in agreement 
with SLHA2 conventions [20]: 

^W = fi'S^ + ^pS, (2.14) 
AVsoft = ml \Sf + [XA^SH^Ha + ^kA^S'' + m'^S' + ^sS + h.c.) . (2.15) 

Then the expressions (12. 7p to (12.121) lead to 

2 

m 



~ M- (2,16) 



(f ~ , (2.17) 



4 

m 



- j^, (2.18) 
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A. = -^A. ~ - (2,19) 
™? ~ ^. (2.20) 

4 

is - — ■ (2.21) 

Next, within the class of models defined by the superpotential (12. ip . there exist the 
diagrams shown in Fig.l which generate terms in AVgoft which are not included in the list 
fl2.16p - fl2.2ip . The corresponding operators and soft terms (after the replacement of Fh^ 
and Fh^ by their tree level expressions) are given by 

2 

^ ^ 9 TTl 

HuH^X + h.c. ^ fr^«^H„ ^ AAa = AA^ ~ — , (2.22) 



M2 



M 

4 



^^H^H: -^Am^-^, (2.23) 



together with analogous expressions with if„ replaced by (and by A^j). 




Similar expressions are also generated by i) the replacement of the shaded bubbles in 
Fig.l by the effective operators (12. 9p and (I2.10p (which generate the soft terms (I2.18P and 
(I2.17P ). and ii) the Renormalization Group (RG) evolution of A\, At, Ah, and from the 
messenger scale M down to the weak (or SUSY) scale Msusy- Whereas this RG evolution 
sums up potentially large logarithms of the form ln(M2/M|[j^y), it does not describe con- 
tributions without such logarithms which serve as boundary conditions for the RG evolution 
at the scale Q'^ = M"^. 

Note that both contributions (I2.22p and (12.230 are generated only at (or beyond) two 
loop order, and are hence suppressed by additional factors additional Yukawa 
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couplings. Compared to the effective SUSY brealcing scale m^/(167r^M), the contribution 
to the A terms (12.221) is negligibly small. However, the contribution (12.231) to Am^ = Am^ 
can be of the same order as the two loop contributions mediated by gauge interactions (see 
appendix A), if A is not too small. Since the contribution (I2.23P to Am"^ = Am^ is typically 
negative, we will subsequently parametrize it in terms of defined as 

Ami = ^< = -^^jS^2^lusY (2.24) 
with MsusY = rn? /M asm. appendix A, and Ah bounded from above by Ah < (Yukawa)^ 

<o{i). 

To summarize this Section, within the class of models defined by the superpotential (12. ip 
one obtains in general, after integrating out the messenger/sequestered sector, an effective 
NMSSM valid at scales below the messenger scale M , which includes 

a) the first two terms in the superpotential (12. ip . 

b) the soft SUSY breaking gaugino, squark, slepton and Higgs masses obtained by gauge 
mediation, which we recall for convenience in appendix A, 

c) additional terms in the superpotential (I2.14p and additional soft terms (12.151) . 

d) additional contributions to the soft SUSY breaking Higgs masses as in (12.241) . 

Note that neither an explicit /i term nor an explicit m| = 5/i term are present at the 
messenger scale M. However, once the above soft terms are used as boundary conditions for 
the RG evolution from M down to Msusyi a term of the form m| H^Hd can be radiatively 
generated in general. (In the appendix B, we recall the /3-functions of the parameters of the 
Higgs sector of the general NMSSM. One finds that a non-vanishing parameter m'g generally 
induces a non- vanishing m\.) 

Depending on the structure of the messenger /sequestered sector, many of the terms 
in (I2.16P - (I2.2ip can be disallowed or suppressed by discrete or approximate continuous 
symmetries. (Exact continuous symmetries forbidding any of these terms would be spon- 
taneously broken in the physical vacuum, giving rise to an unacceptable Goldstone boson.) 
An exception is the term (I2.10p leading to the soft singlet mass term (12.180 . which can never 
be suppressed using symmetries. However, precisely this term is often generated only to 
higher loop order and/or to higher order in an expansion in m/M as expected from naive 
dimensional analysis [4, 18]. Finally we remark that terms of the form SFg + h.c. (which give 
rise to the trilinear soft terms (I2.19P ) will be suppressed if an i?-symmetry is only weakly 
broken in the scalar sector. 

3 Phenomenological analysis 

The purpose of this Section is the phenomenological analysis of various scenarios within the 
class of models defined in Section 2, that differ by the presence/absence of the different terms 
fl236|) to (Km and (K2^ . 

To this end we employ a Fortran routine NMGMSB, that will be made public on the 
NMSSMTools web page [21]. The routine NMGMSB is a suitable generalization of the 
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routine NMSPEC (available on the same web site) towards the general NMSSM with soft 
SUSY breaking terms specified by GMSB, i.e. it allows for a phenomenological analysis of 
the class of models defined in Section 2. It requires the definition of a model in terms of 
the parameter A and the soft SUSY breaking and superpotential terms b) - d) above. Since 
the coupling A at the effective SUSY breaking scale plays an important phenomenological 
role (and in order to allow for comparisons with other versions of the NMSSM as mSUGRA 
inspired), the coupling A on input is defined at an effective SUSY breaking scale Qsusy 
given essentially by the squark masses. The remaining input parameters, notably the soft 
SUSY breaking terms listed in appendix A and in (12.161) - (I2.2ip . are defined at a unique 
messenger scale M. 

The RG equations are then integrated numerically from M down to Qsusy- Additional 
input parameters are, of course, Mz, and also tan/3 (at the scale Mz). Similar to the 
procedure employed in NMSPEC, the minimization equations of the effective Higgs potential 

- including radiative corrections as in [21] - can then be solved for the Yukawa coupling k in 
the superpotential (12.11) . and for the SUSY breaking singlet mass m| (I2.18P or, if m| is fixed 
as input, for ^5. (If specific values for n, m| and at the scale M are desired as input, this 
procedure is somewhat inconvenient. Then, one would have to scan over at least some of the 
other input parameters and select points in parameter space where k, m| or C,s ^ which are 
given at the scale M as output - are close enough to the desired numerical values.) Since the 
gauge and SM Yukawa couplings are defined at the scale Mz, a few iterations are required 
until the desired boundary conditions at Mz and M are simultaneously satisfied. 

After checking theoretical constraints as the absence of deeper minima of the effective 
potential and Landau singularities below M, the routine proceeds with the evaluation of 
the physical Higgs masses and couplings (including radiative corrections as in [21]) and the 
sparticle spectrum including pole mass corrections. Then, phenomenological constraints can 
be checked: 

- Higgs masses, couplings and branching ratios are compared to constraints from LEP, 
including constraints on unconventional Higgs decay modes [22] relevant for the NMSSM; 

- constraints from i?-physics are applied as in [23], and the muon anomalous magnetic 
moment is computed. 

Subsequently we investigate several scenarios, for which many (but different) terms in 
the list ( I2.16P - ( ]2.2ip vanish or are negligibly small. 

3.1 Scenarios with tadpole terms 

The tadpole terms C,f in AIV in (I2.14p and C,s in ^Vsoft in (12.151) will trigger a nonvanishing 
VEV of S. However, as it becomes clear from (12.171) and (12.211) . these tadpole terms - if not 
forbidden by symmetries - tend to be too large: the scale of the soft SUSY breaking gaugino, 
squark, slepton and Higgs masses in GMSB models is given by Msusv ~ /M (together 
with an additional loop factor (IGvr^)"^, see appendix A). Written in terms of M and Msusy-, 
the maximally possible order of magnitude of the supersymmetric and soft SUSY breaking 
tadpole terms are m? MMsusy and ^5 ~ m^/M ~ MM^jj^y- If M > Msusy, 

which will generally be the case, these tend to be larger than the desired orders of magnitude 
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~ Mg^gY and C,s ~ M^jj^y- (This problem is similar to the fi and problem in the 
MSSM with GMSB, see [14].) Hence one has to assume that these terms are suppressed, 
e.g. generated to higher loop order only as in [3], or involve small Yukawa couplings. Let us 
study the latter scenario quantitatively in a simple model [4]: let us assume that the singlet 

superfield couples directly to pairs of messengers 0, (in 5 and 5 representations under 
SU{5)) due to a term 

- 77^0 (3.25) 
in the superpotential W in (12.11) . Then, one loop diagrams generate [4] 



and 



^F = n,^mHn{AyM^) (3.26) 

in agreement with the power counting rules (12.171) and (I2.2ip . (The UV cutoff A appears in 
(I3.26P only if the SUSY breaking Fx is "hard" in the sense discussed in Section 2; otherwise 
the logarithm in (I3.26P should be replaced by a number of 0{1).) 

Below, we consider a mass splitting m? ~ 8 x 10^° GeV^ among the messenger scalars and 
pseudoscalars, and a messenger scale M ~ 10^ GeV. Then, for In [K^ /M'^) ~ 3, a Yukawa 
coupling ?7 ~ 2 X 10"^ generates ~ (150 GeV)^ and ~ —(1 TeV)^. We find that these 
orders of magnitude for and are perfectly consistent with a phenomenologically viable 
Higgs sector. Given the presence of small Yukawa couplings in the Standard Model, and the 
possibility of obtaining additional symmetries in the limit of vanishing r/, we do not consider 
7] ~ 10~^ — 10~^ as particularly unnatural. 

The coupling (13.250 also gives rise to a positive SUSY breaking mass squared 



for the singlet S. Under the assumption of such small values for 77, this term is numerically 
negligible (as well as contributions to Aa, A^, //', A// and two loop contributions to m| 
of 0(mVM2)). 

Hence in the following we will concentrate on models where, among the terms in (I2.16P 
- (I2.2ip and (I2.24p . only and are nonvanishing. (These models are then similar to the 
ones denoted as "nMSSM" in [24]. However, given the present constraints on the soft terms 
we found that a term ~ k in the superpotential (12.11) is required for the stability of the scalar 
potential.) The remaining free parameters are tan/?. A, M, m^/M and ^p'- since m| is fixed 
as input at the scale M (where m| = 0), the equation following from the minimization of 
the potential w.r.t. S can be used to determine ^5. 

Quite generally, there exist two distinct allowed regions in the parameter space, which 
differ how the lightest scalar Higgs mass m^,^ satisfies the LEP bound of ~ 114 GeV: 
a) region A at low tan (3 and large A, where the NMSSM specific contributions to the lightest 
Higgs mass allow for values above above 114 GeV. Low values of tan/5 demand that the 
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messenger scale M is not too large: tan /5 ~ 1 requires ~ at the SUSY scale, but the 
RG equation for differs from the one for by the presence of the top Yukawa coupling 
(which is particularly large for small tau/S). Thus the range of the RG running should not 
be too big, i.e. the scale M should not be too far above the SUSY scale, 
b) region B at large tan /3, where the messenger scale M is quite large (typically ~ 10^'^ GeV) 
resulting in stop masses in the 1.5-2 TeV range. Then the top/stop radiative corrections 
to the lightest Higgs mass can lift it above 114 GeV without the need for NMSSM specific 
contributions. (At large tan /3, A does not increase the lightest Higgs mass; on the contrary, 
large values of A lower its mass through an induced mixing with the singlet-like scalar. Hence, 
A must be relatively small here.) However, in the present context one finds from (13.261) and 
( ]3.27p that such large values for M (with fixed rr? jM ~ 10^ GeV) would require extremely 
small values for 77. For this reason we confine ourselves to region A in the following. 

In region A, the LEP bound on ruh^ requires tan /5 to be smaller than ~ 2, and A larger 
than ~ 0.45. Subsequently we investigate the interval 0.45 < A < 0.6 and tan/? > 1.2, 
where perturbativity in the running Yukawa couplings A, k and ht is guaranteed at least 
up to the messenger scale M. If we naively extrapolate the RGEs beyond the scale M 
(taking the contributions of the messenger fields to the running gauge couplings into account), 
perturbativity in the running Yukawa couplings is usually not satisfied up to the GUT scale 
in region A (in contrast to scenarios where M ~ 10^^ GeV). There exist different possible 
solutions to this problem: first, additional matter could be present at the messenger scale, 
charged under the SM gauge groups. Then, SM gauge couplings can become large (at 
the boundary of perturbativity) below the GUT scale, and since they induce a negative 
contribution to the (5 functions for ht and A, they could help to avoid a Landau singularity 
in the Yukawa sector below Mqut- Another attitude would be to assume that a strongly 
interacting sector (possibly responsible for the breaking of supersymmetry) exists at or above 
the messenger scale M; then the singlet 5, for example, could turn out to be a composite 
state which would imply a compositeness condition equivalent to Landau singularities in the 
Yukawa couplings of S at the corresponding scale (without affecting, at the one loop level, 
the grand unification of the SM gauge couplings). 

Within the region 1.2 < tan/3 < 2 and 0.45 < A < 0.6, a wide range of the remaining 
parameters M, rr? jM and satisfies all phenomenological constraints. Subsequently we fix 
these parameters near the center of the allowed range: M = 10^ GeV, m?/M = 8 x 10^ GeV 
and .^i? = 3 X 10^ GeV^, and vary tan/5 and A in the above intervals (taking, for simplicity, 
ns = 1). 

The allowed range of tan/? (tan/3 < 1.6 for these values for M, m^/M and ^p) and A 
(actually A > 0.5) is shown in Fig. 2; the upper limit on tan/? originates from the LEP 
bound on the lightest Higgs mass mh^. This becomes evident from Fig. 3, where we show 
the range of m^^ (for various values of A, larger values of A corresponding to larger values of 
m/ij) as a function of tan/3. If we would allow for larger values of A (and/or smaller values 
of tan/3), larger values for would be possible. 

In Fig. 4 we display the charged Higgs mass mh± (practically degenerated with a scalar 
with mass m/12 and a pseudoscalar with mass maj), the singlet-like scalar mass m/jg and 
the singlet-like pseudoscalar mass maj, all of which are nearly independent of A. For small 
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0.52 ^■■HiiMkSM'-''' ■ 

0.5 ^ 

0.48 ^ • 

1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 

tan p 

Figure 2: Allowed values of A as a function of tan/3 for M = 10^ GeV, Msusy = Tn?/M = 
8 X 10^ GeV and = 3 x 10^ GeV^. 



tan/5 the large values of the Higgs masses indicate that this region is implicitly more fine 
tuned. The remaining sparticle spectrum is essentially specified by M and m'^/M, and hardly 
sensitive to tan f] and A within the above intervals: 



Bino : 


r\j 


105 GeV 


Winos : 


rsj 


200 GeV 


Higgsinos : 


rsj 


670 - 1000 GeV 


Singlino : 


rsj 


900 - 1800 GeV 


Sleptons : 


rsj 


140 - 290 GeV 


S quarks : 


rsj 


640 - 890 GeV 


Gluino : 


rsj 


660 GeV 



(Due to the small value of tan/5 in this scenario, the supersymmetric contribution to the 
muon anomalous magnetic moment is actually too small to account for the presently observed 
deviation w.r.t. the Standard Model.) 

In Fig. 5, we give the values of C,s (at the scale M), which are obtained as an output 
as function of tan/5. Within the model corresponding to f l3.25p - f l3.27p above, one can 
easily deduce the Yukawa coupling t] from ^5 using f l3.27p resulting in 77 varying in the range 
2 X 10"^ (for tan/5 = 1.6) to 10"^ (for tan/5 = 1.2). The corresponding value of In (A^/M^) 
can then be deduced from (I3.26p . with the conclusion that In (A^/M^) should assume values 
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m. (GeV) 




Figure 3: The lightest Higgs mass as a function of tan (5 for the same parameters as in Fig. 2, 
larger values of rrih^ corresponding to larger values of A. 



in the range 1 to 4 - a reasonable result, by no means guaranteed, that we consider as a 

strong argument in favour of such a simple model. 

Finally we note that for larger values of (as = 3), M (as M = 2 x 10^°) and 
(as = 10^ GeV^, see also the next subsection) phenomenologically viable regions in 

parameter space exist where the running Yukawa couplings A, k and remain perturbative 

up to Mgut- Within the model above, these scenarios would require an even smaller Yukawa 

coupling ?7, T) ~ 10~®. 

3.2 Scenarios without tadpole terms 

Scenarios without tadpole terms have been proposed in [16]. If the number of messengers 

is doubled (77,5 = 2), i.e. introducing $1, $1, $2 and $1, these can couple to S and to the 
spurion X in such a way that a discrete symmetry is left unbroken by the VEV of X [16]: 

W = X + f 2$2) + r/Sf i$2 (3.29) 

Then tadpole terms ~ S,f and ~ ^5 are disallowed, and the Yukawa coupling 77 can be 
much larger. These scenarios have been recently investigated in [18] (see also [19]), where 

the SU{b) breaking (generated via the RG equations between Mqut and M) inside r/S' $i$2 
has been taken into account. 
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For larger values of rj, messenger loops generate non-negligible values for the singlet mass 
m| (12.181) . trilinear A-terms f l2.19p and corrections Am^ = Am^ as in fl2.24p at the scale 
M [18]. Phenomenologically viable regions in parameter space have been found in [18], where 
the parameters M and MgusY have been chosen as M = 10^^ GeV and Msusy = m? /M = 
1.72 X 10^ GeV. The stop masses are quite large (up to ~ 2 TeV) such that the stop/top 
induced radiative corrections to lift it above the LEP bound of ~ 114 GeV. 

We have re-investigated this scenario in a somewhat simpler setup: first we observe that 
the generated values for and Ah, in the notation (I2.19P and (I2.24p . are always related by 

^« = -t^^hMsusy (3.30) 

(with Ah = 2^1) + 3^|. in the notation of [18], where C,d,t denote Yukawa couplings corre- 
sponding to our T] in (13.291) . At Mqut one has = = [18]-) The singlet mass at the 
scale M is then of the order 

ml ^ Qa^, - ^{16gl + 6gl + ^g^An - i^'An^ Mj^sY , (3-31) 

where we have neglected the SU (5) breaking among the Yukawas at the scale M. 

We tried to reproduce the three phenomenologically viable regions in parameter space 
studied in [18]: region I where <^ 1, region III where 0.6 < < 1.1, and region II where 
1-3 ^ 1^(7 ^ 2. We observe, however, that for > 0.7 (or Ah > 1-5 after taking the running 
of between Mqut and M into account) the generated value for \AJ^ from (13.301) exceeds 
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Figure 5: $,s as a function of tan /3 for the same parameters as in Fig. 2. 

~ 5 TeV at M (still > 2 TeV at the weak scale), which we interpret as a certain amount of 
fine tuning between the remaining parameters of the Higgs potential. We will not consider 
the region II below. Note that, as in [18], we obtain Kj clS clll output (from the minimization 
equations of the Higgs potential with Mz as input), which can hide the fine tuning required. 

Limiting ourselves to Ah < 1.5 {\Af^\ < 5 TeV), we were able to confirm the region I. In 
Table 1 we show the Higgs spectrum, and in Table 2 the essential features of the correspond- 
ing sparticle spectrum for a representative point PI in region I, where = —160 GeV, 
Ah = 0.1, A = 0.02 and tan/3 = 6.6 (leading to m| ~ —2.8 x 10^ GeV^ in agreement 
with fl3.3ip ). The point P2 in Tables 1 and 2 is in the region HI of [18]: there one has 

= -4.77 TeV, = 1.46, A = 0.5 and tan/5 = 1.64 (m| 5.3 x 10^ GeV^). We 

see that, in spite of stop masses in the 2 TeV region, rrih^ is not far above the LEP bound. 
On the other hand these results confirm the phenomenological viability of the scenario pro- 
posed in [16,18]. (However, due to the very heavy sparticle spectrum the supersymmetric 
contribution to the muon anomalous magnetic moment is still too small to account for the 
presently observed deviation w.r.t. the Standard Model.) 

3.3 Scenarios without tadpole and 74-terms 

The scenario discussed in the previous subsection belongs to those where many (actually 
most) of the operators (12.71) - (I2.12p and (12.221) - (12.231) are forbidden by a discrete 
symmetry, which is left unbroken in the messenger /sequestered sector, but under which S 
carries a non-vanishing charge. In the above case - where Zjy is not an i?-symmetry - all 
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Table 1: Input parameters and Higgs masses for five specific points. 



Point 


PI 


P2 


P3 P4 P5 


Input parameters 




Messenger scale M (GcV) 


10^^ 


10^^ 


4 X 10** 


3 X 10' 


5 X 10^* 


MsusY = m'/M (GeV) 


1.72 X 10^ 


1.72 X 10^ 


3.2 X 10-* 


3.5 X 10-* 


7.5 X 10* 


tan/3 


6.6 


1.64 


1.6 


1.9 


40 


ris 


2 


2 


2 


2 


2 


A 


0.02 


0.5 


0.6 


0.6 


0.01 




-160 


-4770 













0.1 


1.46 











mi; (GcV") 


-2.8 X 10^ 


-5.3 X 10'' 


-4.3 X lO'' 


-2.1 X lO'' 


-5.0 X 10^ 


CP even Higgs masses 




771^0 (GcV) 


116.1 


115.8 


115.5 


96.1 


94.5 


TO^o (GcV) 


794 


2830 


607 


514 


120 


m-ho (GcV) 


1762 


3411 


717 


579 


603 


CP odd Higgs masses 




m„o (GeV) 


448 


2842 


40.5 


11.5 


1.1 


?7i„o (GeV) 


1761 


3662 


628 


546 


603 


Charged Higgs mass 




m^± (GcV) 


1764 


2862 


619 


535 


613 



Table 2: Some sparticle masses and components for the five specific points of Table 1. The 
chargino masses are close to the wino/higgsino-like neutralino masses, the right-handed/ left- 
handed slepton masses close to the staui/stau2 masses, and the remaining squark masses 
are of the order of the gluino mass. 



Point 


PI 


P2 


P3 


P4 


P5 


Neutralinos 




Xi mass (GeV) 


467 


469 


80.5 


88.3 


101 


Dominant component 


bino 


bino 


bino 


bino 


singhno 


X2 mass (GeV) 


839 


890 


152 


166 


200 


Dominant component 


singiino 


wino 


wino 


wino 


bino 


X3 mass (GeV) 


882 


2322 


463 


428 


380 


Dominant component 


wino 


liiggsino 


liiggsino 


higgsino 


wino 


X4 mass (GeV) 


1432 


2325 


476 


438 


675 


Dominant component 


higgsino 


liiggsino 


higgsino 


higgsino 


higgsino 


X5 mass (GeV) 


1440 


4019 


721 


572 


685 


Dominant component 


higgsino 


singiino 


singiino 


singiino 


higgsino 


Staui mass (GeV) 


692 


693 


100 


103 


260 


Stau2 mass (GeV) 


1100 


1096 


188 


198 


514 


Stopi mass (GeV) 


1931 


1819 


376 


459 


872 


Gluino mass (GcV) 


2389 


2386 


522 


569 


1117 
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soft terms m|, = 3 Ax and the parameter Ah in (12.241) will in general be non- vanishing 
(all others being forbidden). 

The fate of /^-symmetries in the context of gauge mediation has recently been reviewed 
in [25]. In the case of spontaneous breaking within the messenger /sequestered sector [26], 
i?-symmetry violating terms in the effective low energy theory will be suppressed relative to 
i?-symmetry conserving terms. Then, the trilinear terms = 3Ax ( 12.19P will be negligibly 
small. Although the i?-symmetry breaking gaugino masses will typically also be smaller 
than the scalar masses at the messenger scale [25], we will consider in this subsection an 
illustrative scenario which is just a limiting case of the one previously discussed. 

We will investigate the case where the trilinear terms vanish, and where only m| (which 
can never be forbidden by symmetries) assumes natural values at the messenger scale M. For 
simplicity, we will allow for standard gaugino masses (and the usual scalar masses) as given in 
appendix A. Now, the scalar sector of the NMSSM has an exact i?-symmetry at the scale M, 
with identical charges for all superfields. Given that gaugino masses break this i?-symmetry, 
radiative corrections (the RG running between M and the weak scale) induce i?-symmetry 
violating trilinear terms in the scalar sector. If M is not too large or if A, k are small, these 
trilinear terms remain numerically small, and the i?-symmetry in the scalar sector is only 
weakly broken. Given that this approximate i?-symmetry is spontaneously broken at the 
weak scale by the VEVs of H^, and S, a pseudo Goldstone boson (a pseudo i?-axion [27]) 
appears in the spectrum. Light pseudoscalars can lead to a reduction of the LEP constraints 
on m/jj, and have recently been the subject of various investigations [28]. 

In what follows we study the phenomenological viability of such scenarios, which are 
defined by having all terms fl2.16p - (12.211) vanish except for m| (but vanishing A^, A\). 
For simplicity we will also assume that Ah in (12.241) is negligibly small. Then, the model 
is completely specified by A, tan/3 and the scales M and Msusy (recall that k and m| can 
be obtained from the minimization equations in terms of Mz and of the other parameters). 
Again we found that two completely different regions in parameter space are phenomenolog- 
ically viable. 

As before, the first region is characterized by small values of tan (3 (tan (3 <2) and large 
values of A. Relatively large negative values for the soft mass m| for the singlet of the order 
m| ~ —(600 GeV)^ are required at the scale M in order to generate the required VEV 
of S. The mass of the lightest CP-odd scalar varies in the range < < 50 GeV, 
where the larger values are obtained for larger messenger scales M ~ 10^ GeV: then the 
RG evolution generates relatively large values A\ ~ 25 GeV at the weak scale (whereas 
remains very small), and this breaking of the /^-symmetry induces a relatively large mass for 
the pseudo i?-axion. On the other hand, arbitrarily small values for A\ and hence for 
can be obtained without any fine tuning for lower messenger scales M . In all cases we find 
that the lightest CP-even (SM like) scalar hi dominantly decays (with branching ratios of 
~ 80%) into hi flicti, which allows for nih^ < 114 GeV consistent with LEP constraints. 

For given A, m^^ is nearly independent of the scales M and Msusy, but decreases with 
tan/3. In Fig. 6 we show a scatterplot for m^^ as a function of tan/3, which is obtained 
for A = 0.6, varying M in the range 10'' GeV < M < 5 x 10^ GeV and Msusy in the 
range 3.3 x 10^ GeV < Msusy < 4.3 x 10^ GeV. All points displayed satisfy LEP and 
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-B-physics constraints. (We have chosen 77.5 = 2 messenger multiplets, but similar results can 
be obtained - for slightly different ranges of M and Msusy - for 77.5 = 1.) 

(GeV) 




90 : 



1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 

tan p 

Figure 6: as a function of tan/3 for A = 0.6, 10^ GeV < M < 5 x 10^ GeV and 

3.3 X 10^ GeV < Msusy < 4.3 x 10^ GeV. Points where, in addition to all LEP and B- 
physics constraints, the SUSY contribution to the muon anomalous magnetic moment can 
(fails to) account for the presently observed deviation with respect to the Standard Model 
are denoted in blue/darker (gray/lighter) color. 

In the region tan/3 > 1.7 (where ^108 GeV) LEP constraints are satisfied only 
for < 11 GeV, so that ai hh decays are forbidden and the dominant decays of hi 
are hi aiOi ^ 4 r (still requiring irth^ > 88 GeV [22]). For tan/3 < 1.7, the dominant 
decays of hi are hi — > aiOi — > 4 6, in which case LEP constraints allow for ruh-^ as low 
as ~ 108 GeV. The complete theoretically possible range for m^^ is now allowed by LEP. 
(Fixing, e.g., M = 10® GeV, the complete range 1.2 < tan/3 < 1.7 is compatible with 
LEP constraints on the Higgs sector within the above range of Msusy- For smaller tan/3, 
however, the hidden fine tuning becomes quite large.) 

Now, in some regions in parameter space, the supersymmetric contribution to the muon 
anomalous magnetic moment is > 10~^, which accounts for the presently observed deviation 
with respect to the Standard Model. The blue (darker) points in Fig. 6 (which appear only 
for tan/3 > 1.5) satisfy this condition. In Tables 1 and 2 we present the Higgs and sparticle 
spectrum for points P3 (with tan/3 = 1.6) and P4 (with tan/3 = 1.9), which are inside the 
blue region of Fig. 6. 

Another interesting region in parameter space is characterized by large values of tan j3 
(tan/3 > 30) and small values of A (A ~ 10~^), associated with small values of k (k < 10~^). 
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In this case, comparatively small negative values for the soft mass m| for the singlet of the 
order m| ~ —(70 GeV)^ are required to generate the required VEV of S. Due to the small 
values of A and k, Ax and Af^ remain small after the RG evolution from M down to the 
weak scale, leading to a pseudo i?-axion with a mass m^^ < 1 GeV. Now ai is particularly 
light since, for small k, it simultaneously plays the role of a Peccei-Quinn pseudo Goldstone 
boson. However, due to the small value of A, the couplings of ai (with doublet components 
< 10""^) are tiny, and this CP-odd scalar would be very hard to detect; the branching ratios 
hi — s> aiOi are practically vanishing. 

The CP-even Higgs sector is still compatible with LEP constraints if M is very large 
(and MsusY somewhat larger than above), leading to a sparticle spectrum (and At) in the 
1 TeV range such that top/stop induced radiative corrections lift up the CP-even Higgs 
masses. Interestingly, in spite of A ~ 10~^, large values for n^ff = A (S") still generate a 
large singlet /doublet mixing for the two lightest CP-even scalars. As an example, point P5 
(which gives a satisfactory supersymmetric contribution to the muon anomalous magnetic 
moment) is shown in Tables 1 and 2. rrih^ ~ 94 GeV is well below 114 GeV, but the singlet 
component of hi is ~ 88% implying reduced couplings to gauge bosons. The state /i2 with a 
mass m/i2 ~ 120 GeV has still a singlet component ~ 48%. With the help of its nonsinglet 
components, the detection of both states seems feasible at the ILC [29]. Also, the lightest 
neutralino is a nearly pure singlino (with nonsinglet components < 3 x 10^'^), which would 
appear at the end of sparticle decay cascades [30]. 

Throughout this paper we have not addressed the issue of dark matter. Clearly, within 
GMSB models the LSP is the gravitino, but heavy remnants from the messenger sector can 
also contribute to the relic density [31,32]. Its evaluation would require assumptions on the 
messenger /sequestered sector and the reheating temperature after inflation, and is beyond 
the scope of the present work. On the other hand, general considerations can possibly help 
to constrain the large variety of different scenarios found here. 

4 Conclusions 

We have seen in this analysis that the NMSSM can solve the yU-problem in GMSB models in a 
phenomenologically acceptable way. Our starting point was a derivation of the magnitude of 
all possible supersymmetric and soft terms in a generalized NMSSM, that can be radiatively 
generated by integrating out a sequestered/messenger sector with couplings to the singlet 
superfield S. For the phenomenological analysis, we confined ourselves to scenarios where 
most of these terms are negligibly small. Nevertheless we found a large variety of very 
different viable scenarios. 

Scenarios with singlet tadpole terms are acceptable, if the linear terms in 5* (or 5") are 
generated to higher loop order only, or if at least one small Yukawa coupling is involved. 
A simple concrete model [4] with a direct coupling of S to the messengers is viable for a 
Yukawa coupling t] < 10~^. In the case of models with forbidden tadpole terms, as those 
proposed in [16] and analysed in [18], we confirmed the phenomenological viability observed 
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in [18] (at least for the regions in parameter space without uncomfortably large values of 

Quite interesting from the phenomenological point of view are the scenarios with vanish- 
ing A-terms at the messenger scale: these automatically lead to a light CP-odd Higgs scalar 
as studied in [27,28], which plays the role of a pseudo i?-axion. In view of the simplicity 
with which these scenarios can satisfy LEP constraints, it would be very desirable to develop 
concrete models which generate this structure for the effective NMSSM at the scale M. 

Finally we recall that the Fortran routine NMGMSB, that allowed to obtain the results 
above, will be available on the website [21]. With the help of corresponding input and output 
files, further properties of the points PI to P5 as sparticle masses, couplings and branching 
ratios can be obtained. 

Note added 

After the completion of this paper another viable scenario was proposed in [34], in which 
the singlet does not couple to the messenger /sequestered sector, but where the source of 
supersymmetry breaking in the messenger sector is not SU(5) invariant. 
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Appendix A 

In this appendix we summarize the expressions for the gaugino and scalar masses (at the 
scale M), which are generated by gauge mediation under the assumptions that the messenger 
sector involves (5 + 5) representations under SU{b) (additional (10 + 10) representations 
can be taken care of by adding three units to 77,5) with a common SUSY mass M, and F-type 
mass splittings rn? among the scalars and pseudoscalars. The U{1)y coupling ai is defined 
in the SM normalization (not in the GUT normalization). For convenience we define the 
scale MsusY = m? /M and the parameter x = Msusy/M (typically <^ 1). The required one 
loop and two loop functions are [31,33] 

/i(x) = ^((l + x)ln(l + x) + (l-x)ln(l-x)), 
which satisfy /i(x — ^ 0) = /2(x 0) = 1. 
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Then the gaugino masses are given by 

Ml = ^MsusYfi{x)ln5 , 
An S 
do 

M2 = ^MsusYfi{x)n5 , 
M3 = ^MsusYfi{x)n5 , 

and the scalar masses squared by 

The terms '•^^ are present for SU (2) doublets only, and the terms '-^^ for SU (3) triplets only. 
The hyper charges Y are 





UL/di 


Ur 


dR 








Y 


1 

6 


2 
3 


1 

3 


1 

2 


-1 


±1 



Appendix B 

In this appendix we give the one loop /3-functions for the parameters in the Higgs sector of 
the general NMSSM, defined by a superpotential 

and soft terms 

Ko/t = ml\H^\^ + ml\H,\^ + ml\S\^ + {\A^SH.,H,+ '^A^S'' + mlH^H, + m'^S^ + isS 
+htAtQ3H,,T^ - hAhQ^HdB"^ - KA^L^HaL^ + h.c.) , 

under the assumption J2i Yimf = 0, which is always satisfied for GMSB models. 
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^ (2A2 + 2/s:2) 



167r2 

167r2 
1 

167r2 
1 



{6\^Ax + 6k^A,) 



167r2 
1 

(A^Aa + ShlAb + AhlAr + Sg^Mi + 3glM2) 



^ ^X^Ax + Gh^At + /i^A, + ^gfM, + ^g^M^ + ^glM. 

^ ' \2 A , a 1.2 A , 1.2 A , 1.2 A , '''„2,i^ , o „2 n /r , ■^^„2; 



- ^7?M,^ - 3glMl) 
1 



167r2 



(A2(m2 +ml + ml + Al) + 3/i2(m2 + m| + m| + Al) 



+hUml + ml + ml + Al) " f (m^ - ml) - g^M^ - SgjM^ 
^ {2X'iml + ml + ml + Al) + K'{6ml + 2Al)) 

^ (^(6A' + 3/1? + 3hl + hl-gl- 3gl) + 2A/s:m'| 

+l2{2X^Ax - Sh^tA - 3hlAb - hlA, + g^M^ + ^glM^)) 
1 



167r2 
1 



(m'|(2A2 + Ak^) + 2Afi:m^ + Aii\X^Ax + k^A^)) 
{is{X^ + ^^^) + 2A/i(m2 + m^) + 4/€/i'm| + 2Xml{Ax + 2/i') 



167r2 

+2iF{X^Ax + K^A^ + 2fi:m'|(A« + 2/i')) 
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